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Abstract. We consider the convex sets of QO's (quantum operations) and 
POVM's (positive operator valued measures) which are covariant under a gen- 
eral finite-dimensional unitary representation of a group. We derive necessary 
and sufficient conditions for extremality, and give general bounds for ranks of 
the extremal POVM's and QO's. Results are illustrated on the basis of simple 
examples. 



1. Introduction 

The need for miniaturization and the new quantum information technology 1 
has recently motivated a search for new quantum devices with maximum control at 
the quantum level. Among the many problems posed by the new technology there 
is the need of engineering quantum devices which perform specific measurements 
El EJ EJ or particular state transformations — the so-called quantum operations 
O El IE] — which are optimized with respect to some given criterion. In most cases 
such optimal quantum measurements/operations are covariant^ with respect to 
a group of physical transformations. For the case of a quantum measurement, 
" group-covariant" means that there is an action of the group on the probability 
space which maps events into events, in such a way that when the quantum system 
is transformed according to a group transformation, the probability of the given 
event becomes the probability of the transformed event. This situation is very 
natural, and occurs in most practical applications. For example, the heterodyne 
measurement ^| El is covariant under the group of displacements of the complex 
field, which means that if we displace the state of radiation by an additional complex 
averaged field, then the output photo-current will be displaced by the same complex 
quantity. 

In quantum mechanics the probabilities for a given apparatus for all possible 
states are described by positive operator valued measures (POVM)p3, and we will 
say that the measurement is covariant when its POVM is covariant under a unitary 
group representational E]. For quantum operations (QO), on the other hand, 
covariance means that the output of a group-transformed input state is simply the 
transformed output state — a situation again quite common in practice. Typically 
covariance means that the apparatus is required to work equally well on a full set 
of states which is invariant under a group of transformations. For instance, if one 
wants to engineer an eavesdropping apparatus for a BB84 cryptographic scheme 
[141 IT5] that clones equally well all equatorial qubits, then the optimal cloning 
operation must be covariant under the group G = Z4 of ir/2 rotations of the Bloch 
sphere around its polar axis, which is a subgroup of the group of all axial rotations 
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G = U(l) 16 . Similarly, if one wants to engineer a QO which works equally well on 
all pure states, then the operation must be covariant under the full SU(d) group, 
where d is the dimension of the Hilbert space of the quantum system. 

It is easy to see that all POVM's covariant under some group representation 
make a convex set, which describes the complete class of possible covariant appa- 
ratuses. The same obviously holds for group-covariant QO's. Typically in most 
applications the optimization resorts to minimize a concave function on the convex 
set of covariant machines (in quantum estimation theory^) actually such function 
is generally linear), whence the optimal machine will correspond to an extremal 
element of the convex set. For such purpose it is convenient to classify all extremal 
covariant POVM's and QO's, and this is precisely the subject of the present paper. 

For finite dimensional Hilbert space, a characterization of all non-covariant ex- 
tremal QO's was given in Ref. |17j . whereas a characterization of all extremal 
POVM's can be found in Refs. JS| an d El f° r discrete finite probability space. 
On the other hand, no classification of the extremal QO's or POVM's is avail- 
able yet under a covariance constraint, since, as we will see, this constraint makes 
the classification problem much harder. Coincidentally, in many applications the 
optimal QO/POVM is restricted to be rank-one from the special form of the op- 
timization function (this is the case, for example, of optimal phase estimation for 
pure states [21 1201, or of phase covariant optimal cloning of pure states ^1), and 
this has lead to a widespread belief that optimality is synonym of rank-one. How- 
ever, as we will see in this paper, for sufficiently large dimension the extremal 
QO's/POVM's can easily have rank larger than one: this can actually happen for 
optimization with mixed input states, such as in the case of optimal phase estima- 
tion with phase-coherent mixed states |2T|. 

In this paper we provide a classification for finite dimensions of all extremal 
POVM's and QO's that are covariant under a general unitary group representation. 
We will generally consider continuous Lie groups, since then all results will also 
apply to the case of discrete groups as well, with just a little change of notation. 
We provide necessary and sufficient conditions for extremality, along with simple 
necessary conditions, which allow to "sieve" the extremal QO's/POVM's. From 
these conditions general bounds for the rank of the extremal QO's/POVM's easily 
follow as corollaries. 

The paper is organized as follows. In Sect. EJwe briefly review the concept of 
POVM and that of covariant POVM based on the Holevo's theorem(2]. In Section 
EJwe recall the necessary concepts about QO's, including their operator form intro- 
duced in Ref. 22 , which allows to easily classify the covariant QO's as non-negative 
operators in the commutant of a suitable representation of the group. Section 01 is 
entirely devoted to some technical lemmas which will be used in the classification of 
both POVM's and QO's. Finally Sections El and El contains the classification theo- 
rem of extremal group covariant POVM's and QO's, respectively, with some simple 
explicit examples, in particular with application to phase-covariant estimation and 
phase-covariant optimal cloning. 

2. Positive operator valued measures (POVM) 

In the following we will denote by £?(J?T, Jff) the linear space of bounded opera- 
tors from the Hilbert space jff to the Hilbert space J^, and by B(Jf) = B(Jf, Jf) 
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the algebra of bounded operators on JriP . By T\ (Jt?) we will denote the trace-class 
operators on J$?, and by T^~(J$?) its positive elements. 

A general measurement is described by a probability space X equipped with a 
sigma-algebra structure o~(X) of measurable subsets B £ a(X). The measurement 
returns a random outcome x £ X. In quantum mechanics the probability that the 
outcome belongs to a subset B £ a(X) depends on the state p £ T^(J^) of the 
system in a way which is distinctive of the measuring apparatus according to the 
Born rule 



where P is a function on a(X) which is positive-operator valued in B(Jtf), with the 
normalization condition 



Positivity of P is needed for positivity of probabilities for every state p, whereas 
Eq. J5J guarantees normalization of probabilities. In synthesis, P is a positive 
operator valued measure (POVM) on the probability space X. In a sense the POVM 
P represents our knowledge of the measuring apparatus from which we can infer 
information on the state p from probabilities. The linearity of the Born rule Q 
in both arguments p and P is consistent with the intrinsically statistical nature 
of the measurement, in which our partial knowledge of both the system and the 
apparatus reflects in convex structures for both states and POVM's. This means 
that not only states, but also POVM's can be "mixed", namely there are POVM's 
that give probability distributions that are equivalent to choose randomly among 
different apparatuses. 

2.1. Group covariant POVM's. Let's consider now the general scenario in which 
a group of physical transformations G can act on the probability space X. We will 
write gx for the action of the group element g £ G on the point x £ X, and gB 
for the action of g on a whole subset B C X. We will always consider the case in 
which G acts transitively on X, namely for any two points on X there is always a 
group element which connects them. A consequence of transitivity is that X can 
be always regarded as the homogeneous factor space X = G/G x , G x denoting the 
stability group of any point x £ X. 

A POVM P on Jf? for the probability space X is covariant under the unitary 
representation g — > U g of the group G when for every set B £ a(X) one has 



The following general theorem by Holevo 2 classifies all group-covariant POVM's. 

Theorem 1 (Holevo). For square-integrable representations, a POVM P on the 
probability space X is covariant with respect to the unitary representation g — ► U g 
on Ji? of the group G of transformations of X if and only if it admits a density of 
the form 

(4) dP x = U* x EUg x dx, g x £ G : g x x = x, 

where dx is an invariant measure on X, with 2 > in the commutant G' of the 
isotropy group G Xo ofxo, satisfying the constraint 



(1) 



p(B) = Tr[P(B)p] 



(2) 



P{X)=Ije. 



(3) 



UlP(B)U g = Pig-'B). 



(5) 
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with dg invariant measure on G. 

In the case in which the POVM is designed to estimate the group element itself 
g G G corresponding to an unknown transformation U g , then the stability group is 
the identity, whence X — G and the POVM P is covariant if and only if it admits 
a density of the form 

(6) dP g = UlEU g dg, geG 

for any 3 > satisfying the constraint JSJ). The possible seed operators 3 > 
satisfying the constraint © form a convex set. In Section we will classify all 
extremal elements 5 of such convex set. 



3. Quantum operations 

The mathematical structure that describes the most general state change in 
quantum mechanics — such as the evolution of an open system or the state change 
due to a measurement — is the quantum operation (QO) of Kraus Such ab- 

stract theoretical evolution has a precise physical counterpart in its implementations 
as a unitary interaction between the system undergoing the QO and a part of the 
apparatus — the so-called ancilla — which after the interaction is read by means of 
a conventional quantum measurement. We can consider generally different input 
and output Hilbert spaces ^ and J(f , respectively, allowing the treatment of very 
general quantum machines, e. g. of the kind of quantum optimal doners |2.'{l I22j . 
For example in the cloning from one to n copies one has input space Jt? and output 
space Jtf — J{?® n , or its symmetric version Jtf = (Jf?® n ) + for symmetric cloning. 
Within the present paper we will only consider finite dimensional Hilbert spaces. In 
the Heisenberg picture the QO evolves observables, and will be denoted by a map 
^ from B(J^) — > B(Jf). In the Schrodinger picture the QO evolves states, and 
it is given by the dual map : Ti(Jif) — > Ti(Jf), the dualism being determined 
by the equivalence of the two pictures in terms of the trace inner product, namely 
Tr[j?(X)p] = Tr[^£ T (p)X] for all p € 71 (JT) and for all X e B{X). The maps 
^# and are linear completely positive (CP), namely they preserve positivity of 
the input operator for any trivial extension ^ ®/ona larger Hilbert space that 
includes any possible additional quantum system, .y denoting the identity map on 
the additional system. In the Schrodinger picture the CP property physically means 
that the map from 7i to 7i ( J^f) preserves positivity of any input state 
of the quantum system (with Hilbert space -J^) entangled with any possible addi- 
tional quantum system. The map of a QO must also be trace-not-increasing, 
with the trace Tr[^# T (p)] < 1 representing the probability that the transformation 
occurs, and the input and output states being connected as follows 

By denoting with the identity operator on the Hilbert space Jff, we see that 
the trace-not-increasing condition along with positivity of the map are equivalent 
to the constraint 

(8) Jiijx) = Ke B(Jf), 0<K< 1^. 

For finite-dimensional Hilbert spaces it is convenient to represent the maps ^ 
from B{J(f ) — * B(Jif ) as operators Rj( on Jif ® using the following one-to-one 
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correspondence 

(9) R M = Jt T ®J{\l)(l\), J? T {p) = Trj f f[(I je ®p T )R^}, 

where \I) = J2n \ n ) ® \ n ) ls a nxe d vector in J^f ® , {|n) ® |to)} denotes an 
orthonormal basis for Jtif^Jf, and the transposition t for operators is defined with 
respect to the orthonormal basis for B(J^) taken as real. One can easily 

check the correspondence ©, and injectivity follows from linearity. In addition, the 
operator Rj( is non-negative if and only if the map ^# is CP, and the constraint 
(JSJl in terms of the operator K rewrites as follows 

(10) T±x\R*\ = K, 0<K<Ijr. 

The positive operators Rj( satisfying the constraint Hl()(l make a convex set, which 
is the operator counterpart of the convex set of the corresponding QO's . 

3.1. Group covariant CP-maps. We call the map ^# from B(J^f) to £>(Jf ) 
G-covariant, when 

(11) JK{V^XV a ) = U^(X)U g , Vg e G, 

{U g } and {V^} denoting unitary representations of G over the input and output 
spaces Jif and Jff, respectively. The Schrodinger picture version of identity is 

(12) ..# T {U aP Ul) = V a jr(p)Vj, V 5 eG, 

where Jl T goes from Ti(J^) to Ti (Jf). 

The operator form R^ for maps simplifies the classification of QO's that 
are covariant under a group G, resorting to the Wedderburn's decomposition of the 
commutant of the representation. It is easy to show that the map ^# is G-covariant 
(i. e. it satisfies Eq. if and only if its corresponding operator Rj( is invariant 

under the representation V g ® J7*|22]- In fact, from Eq. 10 using invariance of 
partial trace under cyclic permutation of operators acting only on the traced space 
one has 

=^ T (p) - V^^ T {U g pUl)V a 
' " ' = Tr^{(/^ ® p T )[R^ - (Vj g, U g T )RMV g ® U* g )}}, 

and, since Eq. is a one-to-one correspondence between maps and operators, one 
concludes that 

(14) [R^, V g ® U* g ) = 0, V.geG. 

Therefore, the problem of classifying covariant CP-maps resorts to that of classify- 
ing positive elements of the commutant of the representation V g ® U * on <W (g> ffl 1 . 
By labeling with k the generic equivalence class of the representation, with multi- 
plicity mfe, the Wedderburn's decomposition of the representation space is written 
as follows [24] 

(is) jt® je = ®(jffc<g>c rofe ). 

k 

Then, since Rj( must be a positive operator in the commutant of the representation 
it must have the general form 

(16) Rm = ® k (Ije k ® w\w k ) = W^W, W = ® k (Ije k ® w k ), 
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where u>k is any operator on C mfc , i. e. a nik X m& matrix. Therefore, the classifica- 
tion of covariant trace- not-increasing QO's with ^#(7j^) = K < Ijj? is equivalent 
to classify the operators Rj( of the form (|16|l with the constraint 

(17) Tr ^[( J * ® <4™fc)] = K < J>. 
fc 

The constraint (|17|l is generally quite involved, due to the subspace mismatch be- 
tween the tensor product J(f <g> Jtf and the Wedderburn's decomposition: its sim- 
plification will be the main task of Section [BJ 

4. Technical lemmas 

This section will be entirely devoted to technical lemmas, which will be used 
for the classification of both extremal covariant POVM's and QO's. The lemmas 
connect conditions on the vanishing of partial traces with linear spannings. 

In the following we will make use of the following simple fact for any linear 
space Jzf and a subspace C Jzf : if the only vector of Jzf that is orthogonal to 
the whole subspace 5? is the null vector, then one has 5? = Jzf. Moreover, since 
orthogonality to a set s of vector implies orthogonality to its linear span Span(s), 
then the previous assertion holds also for subsets s C Jzf (not necessarily subspace), 
namely if the only vector orthogonal to the subset s is the null vector, than one has 
Jzf = Span(s). From now we will also make use of the following natural notation 

(18) X(B(s*) <8> I a )Y* = Span{X(A <g> Iss)Y\A e B(*f)}, 
for X, Y any operators with domain stf ® 3§. 

Lemma 1. Let B £ B(S§2 <8> srf and 38\ t i denoting arbitrary finite di- 

mensional Hilbert spaces. Then, the injectivity of the linear CP map W{A) — 
Tiag 1 [B^ AB] on B{s?f) is equivalent to the spanning condition 

(19) £(0 = B{B{Sh) ® I®i)B^. 

Proof. The iniectivity of the map W{A) = Tr S g 1 [B^AB] on B(sf) means that 

(20) VA e B(s/) Tisg, [B^AB] = A = Q. 

The condition Tr^B^AB] = is equivalent to Tr[C Tr^ 1 [B^ AB}} = VC € 
B{S§2)- Therefore, since one has 

(21) Tr [CTr^ 1 [5^ AB] ] = Tr[(C <g> I ai )B*AB] = Tr[B(C <g> /.^ )B t A] 
condition (|20|l is then equivalent to 

(22) VA e B(jzOi Tr[B(S(^ 2 ) ® /*)-BU] =0 =^ A = 0, 

where we used notation (|18p. Eq. (|22|l says that the only operator A G B{srf) 
orthogonal to the operator space B{B{S§2) ® Igg^Bt C S(,e/) is the null operator, 
which means that B{JB{SS-i)® Iss-^B^ is actually the full linear space B(£/), namely 
condition (|22|l is equivalent to condition l|19|) • B 

The above theorem leads immediately to the following corollaries. 

Corollary 1. A necessary condition for injectivity of the map W{A) = Tr$g 1 [B^ AB] 
on B{stf) is 

(23) dim(j^) < min{dim(^ 2 ),rank(S)}. 
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Corollary 2. The injectivity of the map W(A) — Trgg^Bi AB] on B(s?) is equiv- 
alent to the existence of a linear injective map V from B(stf) to B(SS) such that 

(24) VA G B{sf) B(Y(A) <8> I ai )5 f = A. 
The relation between the maps W and 'f is given by 

(25) W{A) = T r3dl [B^B(y(A) ® I m )B^B]. 

Proof. The spanning condition (|19fl — equivalent to the injectivity of the map 
W(A) = Tr S g 1 [B^AB] on Z?^)— guarantees that for each A G B(srf) there exists 
an element, say Va, of B(8@) such that B(Va ® Igg x )B* = A. Consider now an 
orthonormal basis Aj for B(gtf), and denote by Vj any element of of B(3$) such that 
B(Vj®IsgjB^ = Aj. It is clear that the {Vj} can be chosen as linearly independent. 
Now, for every element A e define ^(A) = ^ Tr[AtA]^-. This map is 

clearly linear and injective. The map i / (A) corresponds to a nonorthogonal change 
of basis (from {Aj} to {Vj}) which compensates the nonorthogonal change of basis 
B(Vj O i^J-Bt = Aj, Eq. follows by substituting Eq. JjE) into the map ■ 

We have also the additional lemma. 

Lemma 2. As in Lemma^\ the injectivity of the map W(A) — Tr ^[B'AB] on 
B(sz/) is equivalent to the linear independence of the set of operators {WjWj}, where 
Wi G B(3S\ 1 ^2) are defined from the singular value decomposition B — |V*)(Wi| 
through the identity \Wi) = (Wi 8> i^)]/), \I) G -S§f 2 denoting the fixed vector 
I J) = |Z) /or {|Z) (g> |m)} arbitrary orthonormal basis of 3$f . 

Proof. 

First, notice that the identity \X) = (X (8 7^)17) sets a bijection between 
vectors |X) € <8 38\ and operators X G B(3S\^B§2)- Then, using the singular 
value decomposition B = £\ |Vi)(Wj|, with |V^) G and |Wj) G ^ 2 ® ^1, the 
partial trace in Eq. I|20fl becomes 

(26) Tt mi [BtAB\ = 53(V5|A|^) Tr^ [\Wi){Wj\] = ^(^|A|^> 

where r denotes the transposition for which (X ® = [Jss 1 ® X T )\I), and * 

denotes complex conjugation, i. e. X* = (X T )* . By taking the complex conjugate 
of the last equation and introducing the matrix A,j = (Vi|A|Vj)* G Mat(C) where 
N = rank(B) (iV 2 is the cardinality of the set {W^Wj}), the statement (EOJ is 
equivalent to 

(27) {Aij} G Mjv(C), AijWjWj = =>■ A 4j = 0, Vi, j, 

namely the operators {W}Wj} are linearly independent. ■ 

In the following we will need the following generalization of Lemma ^ 

Lemma 3. Let B G B(®k(^2 ® SS^),^), and denote by Pk the orthogonal 
projector over SB^ ® SS\ , srf and 3s^\ being arbitrary finite dimensional Hilbert 
spaces. 
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The following implication 

(28) A e Tr^w [P k B^ABP k ] = Vfc =^ A = 0. 
is equivalent to 

(29) B{^)=Spzn{B[® k {B{^ ) )®I^ ) )]B^}, 
and necessary conditions are 

(30) dim«) 2 < ^dim(^ fe) ) 2 , 

(31) dim(^) < rank(P). 

Proof. The condition Tr [PfcP^APPfc] = Vfc is equivalent to say that for any 

a*! 

C fc € one has Tr[P k C k Tr^w [PtPUPPfc]] = Vfc. Since one has 

Tr[C fc Tr w [P fe pt ABP k \] = Tr[(C fc ® 7 w )P k B^ABP k ] 

(32) 

= Tr[5P fc (C fe ®7^))P fe 5U], 
and, therefore, condition (|28|) is equivalent to 

(33) A e B(jjO, Tr[BP k (B(^ { 2 k) ) <g> 7^, )P fc P t A] = OVfc =► A = 0. 

The last condition says that the only operator in B(jrf) which is orthogonal to the 

set BP k (B(3§^) <g> Iag(k))P k B^ Vfc is the null operator, or, in other words that the 

set spans the full operator space B(£f), namely Eq. 129|) . The necessary conditions 
then follow trivially. ■ 

We are now ready to classify the extremal group covariant POVM's and QO's in 
the following sections. In order to classify extremal elements of convex sets, we will 
use the method of perturbations. We will call a non null operator B a perturbation 
for an operator A in a convex set if both A±tB are still in the convex set for some 
(sufficiently small) t > 0. Then, clearly A is not extremal in the convex set if and 
only if it has a perturbation. 

5. Extremal covariant POVM's 

We have seen that the covariant POVM for the estimation of a group element g 
of an unknown unitary transformation U g is of the general form 

(34) dP g = dgU$EUl 
with probability space X = G, and with 

(35) / dgUpU g =I^. 

Jg 

The Wedderburn's decomposition Ijl5|l of the representation space here rewrites as 
follows 

(36) Jf = 0(jr fc ®C mfc ), 

fe 

where we remind that fc labels the equivalence class of irreducible components, 
and m k denotes its multiplicity. The integral in the normalization condition l|35|l 
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belongs to the commutant of the representation, whence it can be rewritten as 
follows 



(37) / dgUlEU g = ^dJ k [l^ k ®TT^ k (P k EP, 

Jg k 



P k denoting the orthogonal projector on the subspace Mk <8> C mfe . Eq. (|37|l follows 
from the simple fact that for an irreducible representation on the space say Jzf, one 
has J G dgW g ZU g = d^TT[Z]I^ for measure dg normalized to unit on G. Eq. 
(|37|l allows to split the constraint (|35|l into the following set of constraints 

(38) Tr^ fe (P fe SP fc ) = d^J mk , Vfc, 

where by I mk we denote the identity matrix over C mk . We then conclude that the 
classification of extremal G-covariant POVM's is equivalent to find the extremal 3 
within the convex set of operators 3 > satisfying the constraints l)38[l. For such 
purpose we have the following theorem. 

Theorem 2. Let 3 be an element of the convex set of positive operators on ,J%? 
satisfying the constraints 

(39) Tr^ fc {P k EP k ) = d^ k I mk , Vfc G S, 

where S denotes the set of equivalence classes of irreducible components in the rep- 
resentation. Write 3 in the form 3 = X'AX with A > 0, choosing Rng(X) = 
Supp(A) = Ker(A)- L . Then 

1 . O is a perturbation of 3 if and only if is Hermitian, with Tr ^ (P k OP k ) = 
Vfc G S, and = X^BX for some nonzero Hermitian B with Supp(_B) C Supp(A). 

2. For the specific choice of the form of A as A — ® k A k , with A k G B(Jf k <&C mk ), 
one has B = ® k B k , B k G B{J*t? k <g> C mfc ) and Supp(P fe ) C Supp(Afc), Vfc G S; 

3. 3 = X^X is extremal if and only if 

(40) B(RngpQ) = Span{X[® fe (/^ (fe) ® B(C m "))]X^}. 



Proof. 

1. Let 6 Hermitian, with Tr^ fc (P fe 0Pfc) = 0, and 9 = X^BX for some nonzero 
Hermitian B G B(Jf ) and with Supp(P) C Supp(A). Then for rank(P) > 
is necessarily nonzero, and since A > 0, both constraints A ± tB > and 
Tr i ^ fc (Pfe(5 ± tQ)P k ) = d^ k l mk Vfc are satisfied for some t > 0, whence is a 
perturbation for 5. Conversely, suppose 9 G B(J?°) is a perturbation for 3. Since 
we must have 3 ± t@ > and Tr^ fc [P k (E ± tQ)P k ] — d^ k I„ lk for some t > 0, then 
is Hermitian with Tr j%> k (P k QP k ) = Vfc G S. Moreover, if we write 3 in the form 
5 = X^ AX with nonnegative A G B(J^f), and Rng(X) = Supp(A), then also can 
be written in the same form = X^ BX for some nonzero Hermitian B G B(J^) 
and Trj? k [-Pfc (3 ± tO)P k ] = d^ k I mk . In fact, if X is not invertible, it can be always 
completed to an invertible operator Z = X + Y by adding an operator Y with 
Rng(Y") = Ker(A), and one can equivalently write 3 = Z'AZ. Now we can write 
also the perturbation operator in the form = Z^ BZ. However, since A ± tB > 
for some t, then necessarily B must have Supp(P) C Supp(v4) = Rng(X), whence 
Z^BZ = X^BX. 

2. First it is obvious that a choice of the form A = ® k A k , with A k G B(Jif k ® 
C" lfc ) is always possible. Then, in order to have A ± tB > for some £ > 0, one 
must have B = (B k B k , each B k Hermitian, with Supp(Bfc) C Supp(Afc), Vfc G S. 
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3. Since Supp(A) C Rng(X) and A > 0, we can always merge \f~A into X by 
substituting X — > \ AX. Then, since S is not extremal iff it has a perturbation, 
by part 1 one sees that S is extremal iff for Hermitian B € B(Jtf) with Supp(i?) C 
Rng(X), one has 

(41) Trje k (P k X^BXP k ) = 0\fkeS B = 0, 

whence via Cartesian decomposition of B we have the equivalent statement 

(42) B e S(Rngpf)), Ti^ k {P k X^ BXP k ) = OVA; e S => 5 = 0. 
Then, by Lemma lathis is equivalent to condition Q40JI. ■ 

Corollary 3. ^4 necessary condition for extremality of the seed 2 of a group co- 
variant representation as in Theorem\^is 

(43) rank(~) 2 <^m 2 . 

k 

Proof. Eq. 143(1 is a trivial consequence of the necessary condition 1(40(1 . ■ 
Corollary 4. Every rank-one POVM is extremal. 

Proof. For rank(V) = 1 the iff condition 1(40(1 is trivially satisfied. ■ 

Theorem 3. For S containing only a single equivalence class, say h, with multiplic- 
ity nih > 1, the extremality of a covariant POVM on the Hilbert space M' = ® 
C mh is equivalent to the linear independence of the set of operators {W^Wj}, where 
Wi G B(C nih ,J#h) are defined from the spectral decomposition 5 = YV \W)(Wi\ of 
the seed 5 of the POVM through the identity \Wi) = (W, ® I mh )\I) , \I) € (C m ")® 2 
denoting the fixed vector \I) = |Z) ® \l) , for ® \m}} arbitrary orthonormal 
basis of (C" 1 ' 1 )® 2 . Extremal POVM's with any rank rank(S) < m/j are admissible. 

Proof. For S containing a single equivalence class h with multiplicity m/, > 1 the 
seed S of the POVM must satisfy the single constraint 

(44) Tr^ i (S)=^7 mh . 

Now, write S in the form H = X^AX with X € B{JT h <g> C mh ,srf), and Rng(V) = 
Supp(A), &/ being a Hilbert space such that Supp(^4) C srf C ^ l (g)C m '', and which 
can be chosen as ~ Rng(V). Then, according to Theorem[21© is a perturbation 
for 5 iff it is of the form O = VTBV, with £? Hermitian, Supp(£>) C Supp(yl), 
and Trjjf h (X> BX) = 0. This means that the extremality of S is equivalent to the 
injectivity of the map W(B) = Tr^X^BX) over the set of Hermitian operators 
B with Supp(-B) C Supp(A), which is equivalent to injectivity of the same map 
on B(Rng(X)). We are thus in the situation of Lemma with = Rng(V), 
S3\ = C mh and S$i = J^. Therefore, by writing the singular value decomposition 
of X = J2i \Vi)(Wi\, with Span{|Fi)} = Rng(V) = Supp(^l) the injectivity of the 
map W(B) — Tr ^ h \X^BX\ on $(Rng(i2T)) is equivalent to the linear independence 
of the set of operators {W}Wj}, where Wi € B{C mh ^ h ) are defined through the 
identity \Wi) = {Wi ® I mh )\I), \I) G (C m ")® 2 denoting the fixed vector \I) = 
J2l \l) ® with {\l) ® |m)} arbitrary orthonormal basis of (C" 1 ' 1 ) 8 ' 2 . Now, the 
maximum rank of the POVM is given by the maximum number of operators Wi 
such that the set of operators {W}Wj} in B(C mh ) is linearly independent. Since 
we can have at most m 2 linearly independent operators in B(C mh ), the maximum 
cardinality of the set {Wi} is rut - ■ 
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Corollary 5. A POVM which is covariant under an irreducible representation is 
extremal iff it is rank one. 

Proof. For S containing a single equivalence class h with multiplicity rrih = 1 the 
iff condition Q40(l rewrites 

(45) #(Rngpf )) = Span{X(/^. w ® C 1 )^} = Span{ W 1 }, 

which is satisfied iff rank(V) = 1. As an alternative proof, the present corollary 
corresponds to the situation of Theorem |3] for multiplicity = 1. ■ 

5.1. Example. Consider a POVM on Jrff with dim(Jff) = d covariant under G = 
U(l), with 

d-l 

(46) Lfy = exp(i(/)N), N = ^n|ra)(ra|. 

n=0 

Here we have d one-dimensional irreducible representations with characters \k {4 1 ) = 
exp(ik<p), k = 0, ...d — 1, namely they are all inequivalent, whence with unit 
multiplicity. Therefore, the necessary condition (|43|) bounds the rank of the POVM 
as follows 

(47) rank(S) 2 < dim(Jf), 

and in order to have rank(S) = 2 one must have dim(Jff) > 4. According to 
Theorem |21 the extremal POVM's have seed of the form 5 = X'X satisfying the 
identity 

(48) S(RngpO) = Span{\X k )(X k \ : < k < dim(Jf)}. 

where \X k ) — X\k), {\k)} denoting any orthonormal basis for Jf?. Notice that in 
the present example the operator S corresponds to a so-called correlation matrix, 
namely a positive matrix with all ones on the diagonal. This follows from the con- 
straint H38|) . which in our case is simply (fc|S|fc) = 1, Vfc. Therefore, the present 
classification of extremal POVM's coincides with the classification of extremal cor- 
relation matrices given in Ref. |25j . 

5.2. Example. Consider a POVM for n qubits on the Hilbert space Jf = (C 2 )®" 
covariant under the tensor representation U? n of G = U(l), with 

(49) ^ = exp(#|l)<l|), 

where {|0), |1)} is a orthonormal basis for C 2 . Here we have n + 1 one-dimensional 
irreducible representations with characters Xk(4>) = exp(ik(p), k = 0, . . .n, and with 
multiplicity mj, = (Z). An orthonormal basis of each subspace C mfc of — (Bk < C mk 
is given by 

(50) {\j) k } = { P ^\ oo^m^i)}, 

n—k k 

where pj- n,fc - ) denotes the jth permutation of k qubits in the state |1) in the tensor 
product of n qubits in total, with all other qubits in the state |0). In the present 
example, the iff condition for extremality (|40|l requires that 5 = X^X satisfies the 
identity 

(51) B(RngpQ) - Spzn{X\i) kk (j\X\ k G S, i,j = 1, . . . m fe }, 
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where now {|j)fc} denotes any orthonormal basis for C mfc . The necessary condition 
(gSJl bounds the rank of the POVM as follows 

, 52 , ^<-±$=( 2 :). 

Here, in order to have rank(S) > 2 one needs n > 2 qubits. For n = 2 according to 
the previous example, one necessarily must have at least two inequivalent classes, 
since each of the irreducible components has less than four dimensions (the same 
is true also for n = 3). The previous example is also recovered by considering 
the special case in which Rng(X) C ((C 2 )®")+ i. e. containing only the sub- 
representation of U? n on the symmetric subspace ((C 2 ) lg, ™)+, with multiplicity 1. 

5.3. Example. Consider a POVM on ,34?® 2 which is covariant under the group 
representation U a ® Ijg>, where U g is an irreducible representation of G on Jf? . 
Here, we trivially have a single equivalence class, say h, (corresponding to the 
irreducible representation U g ) with multiplicity — dim(^), i. e. the Hilbcrt 
space J$? coincides with the multiplicity space Jf? ~ C mfl . This is exactly the case 
considered in Theorem Therefore, the extremality of the POVM is equivalent 
to the linear independence of the set of operators {W^Wj}, where Wi G B(Jf) 
are defined from the spectral decomposition S = TV |Wi)(Wj| of the seed S of the 
POVM through the identity \Wi) = {Wi <g) as in Theorem Therefore, 

we can have extremal POVM's with any rank(S) < dim{Jf). Notice that there 
cannot be more than a single maximally entangled vector | Wi) in the decomposition 
of S, since, otherwise, at least two operators Wi would be proportional to unitary 
operators, and then the set {w/ W,-} would be necessarily linearly dependent (two 
products would be both proportional to the identity). The rank-one case with a 
single maximally entangled projector corresponds to a so-called Bell POVM. 

6. Extremal covariant quantum operations 

In the following we will denote shortly by .Ag the operator algebra generated 
by the group representation V g (g> U*, by A' G its commutant, and finally by Ti.' G 
the Hermitian operators in the commutant. The following theorem classifies all 
extremal G-covariant maps ^# in the convex set given by Eq. (|17|l . 

Theorem 4. Let R be an element of the convex set of positive operators in the 
commutant A! G of the operator algebra Ac generated by the group representation 
V g <g> U* onX ® Jtf, i. e. of the form 

(53) R = ®k(Ij? k <8 wlw k ) = W^W, W = ® k (Ije k ® w k ), 
satisfying the constraint 

(54) Tr ^[( J « ® w l w k)} =K< 

k 

where 

(55) ® x = [M ® c mk ) 

k 

is the Wedderburn's decomposition of the representation space, k labeling the equiv- 
alence class of representations, with multiplicity nik- Denote by Pk the orthogonal 
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projector over the space Mk <8> C mk of the equivalence class. Write R in the form 
R = X^QX with Q,X G A' G and Rng(X) = Supp(Q). Then: 

1. S is a perturbation of R if and only if S £ Hq, with Tr^[5] = 0, and 
S = X^OX for some nonzero O E Ti.' G with Supp(O) C Rng(X). Specifically, 
writing Q = ®k(Ij? k ® Qk) and X = ®k{I.j^ k ® X k ), one has O = ©fc(Ij^ ® O k ) 
with Supp(O fe ) C Rng(X fe ) Vfc. 

2. One can always write R in the form R = X'X, with X G A' G of the form X = 
®k{Ijv k ®Xk). Denote by S the set of equivalence classes k for which Xk ^ 0. Then, 
a necessary and sufficient condition for extremality of R = X^X with Trjg-[R] = K 
is the injectivity of the map S^iO) — Trj^[X^OX] on A' G (1 B{Rng{X)) , namely 

(56) O G A' G H B(Rng(X)), Tijt[X ] OX] =0=^0 = 0, 
which is equivalent to 

(57) ® feeS S(Rng(X fc )) = (B ke sX k Tr^ [P h {I x ® B{^))P k \x\. 
Proof. 

1. Let 5* G Ti. G , with Tr^S 1 ] = 0, and S = X^OX for some nonzero Hermitian 
O with Supp(O) C Supp(Q). Then for rank(O) > S £ H' G is necessarily nonzero, 
and since H' G 3 Q > 0, all constraints: Q±tO G H' G , Q±tO>0, and Tr*- [R ± 
tS] = K are satisfied for some t > 0, whence S is a perturbation for R. Conversely, 
suppose that S € ® is a perturbation for R. Since we must have TL' G 3 
R ± tS > and Tr,^ [i? ± tS] = K for some t > 0, then Se^ with Tr jr[S] = 0. 
Moreover, if we write R in the form R = X'QX with Rng(X) = Supp(Q), then also 
S can be written in the form S = X^OX for some nonzero Hermitian O € W' G . 
In fact, if X is not invertible, it can be always completed to an invertible operator 
Z = X + Y by adding an operator Y G A' G of the form Y = (Bk(I^f k ® ife) with 
Rng(Yfc) = Ker(Qfc) (where Q = ©feCkes, ® Qfc))i an d one can equivalently write 
R = Z^QZ with Q G and Z G -4' G . Now we can write also the perturbation 
operator in the form S = Z^OZ. However, since for some t the operator Q±tO > 
must belong to the commutant A' G , then necessarily O G H' G and Supp(O) C 
Supp(Q) = Rng(X), with Z^OZ = X^OX. Specifically, writing Q = ® k {I^ k ®Qk), 
one has O = ®k{Ijtf k ® O k ) with Supp(O fc ) C Supp(Q fe ) = Rng{X k ) Vfc. 

2. As in part 1 we can always take Q as the identity, and redefine X — > y[QX, 
since Q > 0, keeping X of the form X = ffifc(/j>r fc <8> Jffc), since both operators in the 
product VQ-^ belong to the algebra A' G . From part 1 we then see that R = X^X 
with X G _4' G is extremal if and only if 

(58) O G W G n S(Rng(X)), ^[X^X] = =^> O = 0, 
and via Cartesian decomposition this is equivalent to 

(59) O e A' G n B(Rng(X)), Tr^^OX] =0=^0 = 0. 

Since O € ^4' G n B(Rng(X)) can be decomposed as O = ®k{I.^ k ® Ok) with Ok € 
S(Rng(Xfe)) Vfc G S, then the statement (|59|l is equivalent to 

Vfc g S O k g B(Rng(Jf fc )), 

(60) ^ Ttjr[(Ijft <8 ^fe) f (^ ® O k )(Ijt? k ® = => O fc = OVfc G S, 
fees 
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or else 
(61) 

Vfc e S O k e B(Rng(X fe )), 

Trjt[®kes(Ijsf k O X k )\l^ k ® O fe )(ijr fc <8> X*)] = 0=^O k = OVfc e S, 

The vanishing of the partial trace can be written as the vanishing of the trace 
Tr[® keS (I^ k ® X k )\l^ k (g> O k )(Ijff k <g> X fc )(7jr $ C)] for any C € S(JT), namely 
the vanishing of Tr{(£ k esO k X k Tr^ [ft(Jjr <8 C)P k ]x\} for any C £ and 
upon defining S = ©fcesOfc, the statement (|61|1 rewrites 

S e 8fc eS S(Rng(X fe )), 

(62) , 

Tr{5 ® ke s X k Trjgj, \P k (l x ® B{J?))P k ]xl} = =>■ 5 = 0, 

namely, since the only operator in the linear space (B k esB(Rr\g(X k )) orthogonal to 
the subspace Q> k€ sX k Tr^ k \P k {Jx ® B(Jf))P k ]Xl is the null operator, one has 

(63) ® keS B{Rng{X k )) = ® ke sX k Tv^ k [P k (I^ ® B(^))P fe ]X^. 
■ 

Corollary 6. As in Theorem^ a necessary condition for extremality is 

(64) ^rank(X fc ) 2 < dim(Jf) 2 , 

fees 

Corollary 7. ^4ny rank-one covariant QO is extremal. 

Proof. For rank(X) = 1 the set S must contain only one equivalence class, and 
the iff condition (|57|l of Theorem 01 is then trivially satisfied. ■ 

Corollary 8. For an irreducible representation any extremal covariant QO must 
be rank-one. 

Corollary 9 (Choi). In the non covariant case, a QO j$ from B(Jf) to B(Jf) 
is extremal iff it can be written in the form ^#((9) = 5Z f Wj OWi, with Wi £ 
B(J$f,Jff) and the set of operators {Wj\Vj} linearly independent. 

Proof. The non covariant case corresponds to the trivial covariance group G = I, 
i. e. the group containing only the identity element. This corresponds to have just 
a single equivalence class, with multiplicity equal to dim(Jf ® Then, as in 

the proof of point 2. of Theorem H the extremality of R = X^X £ B(Jf <g> JT) 
is equivalent to the injectivity of the map W (A) = Tr,j^[X' AX] on B(Rng(X)). 
According to Lcmma|2 using the singular value decomposition X = \Vi){Wi\, 
with \Vf) orthonormal basis for Rng(X) and \Wt) £ ® <ffl , one has ^#(0) = 
W}OWi for O £ B(Jt), and W{A) = Y.^WiWlW* for A £ B(Rng(X)), 
and injectivity of W is equivalent to linear independence of the set of operators 

{w}w 3 }.m ' 

Corollary [5] is the same as Choi theorem ^7]. Notice that differently from the 
case of QO's, for POVM's the non covariant case cannot be recovered as a special 
case of the covariant classification, since the group itself (or, more generally, the 
homogeneous factor space) coincides with the probability space X of the POVM, 
whence trivializing G also trivializes X. 
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6.1. Example. Consider the phase-covariant eloning [TT)l I22| for equatorial qubits 
from 1 to 2 copies. This correspond to G = U(l), with representations U<j, — 
e ^|i)<i|o an d y^ = e #Ea=i where s = denotes the input qubit and s = 1, 2 

the output ones. Here J4? = C 2 and = J^" 8 " 2 . We first need to decompose 
the representation V$ (8> J7f. This is made of one-dimensional representations, 
with characters e lfc< ^, with fe = —1,0,1,2 and multiplicities m_i = 1, mo = 3, 
mi = 3, and rri2 = 1. The necessary condition 164J1 in the present case becomes 
SfceS rank(Xfc) 2 < dim(.Jf) 2 = 4, which means that we can have either a single 
equivalence class with rank(Xfc) < 2, or two equivalence classes with rank(Xfc) = 1 
each. Orthonormal bases for the supporting spaces Jf?k (g> C mfc = C mfe of the fcth 
equivalence class of irreducible representations are reported in Table ^ as subset of 
an orthonormal basis for the tensor product Jtf <g> Jt? '. 



k 


\h) ® |^-) 


-1 


|001) 





|101), |011), |000) 


1 


|100), |010), |111) 


2 


|110) 



Table 1. Orthonormal bases for the supporting spaces J$?k <g> 
C" lfc = C" lfc of the fcth equivalence class of irreducible representa- 
tions for 1 to 2 phase-covariant cloning. The orthonormal basis are 
chosen as subsets of an orthonormal basis for the tensor product 

jt ® jr. 



s = W 


K>>} 






{-1,2} 


1001) 


|no) 




{0,1} 


a|000>+6|011> + c|101) 


a'|lll) +6'|100) +c'|010) 


|a|^ + |6'| 2 + |c'| 2 = l 
|a'| 2 + |£*| 2 + |c| 2 = l 


{o,-i} 


|000> + o]011> +6|101) 


c]001) 


|a| 2 + |6| 2 + |c| 2 = l 


{1,-1} 


a|100) + &|010) + cjlll) 


djOOl) 


|a| 2 + |6| 2 = l 
|c| 2 + |d| 2 = l 


{1,2} 


a|100) + 6|010) + | 111) 


d]110) 


|4* + |6| a + |d| a = l 


{0,2} 


a|000> + 6|011> + cjlOl) 


d]110) 


\a\ l + W = l 
|6| 2 + |c| 2 = l 


{0} 


_^|101) + -^| 11),|000) 






{1} 


_^|oio) + -^1100),! in) 







Table 2. Cloning from 1 to 2 copies: classification of operators 
R = J2kes ^ k = lW ) I satisfying the necessary condition. 



The operators R = J2keS^ k = m /W I satisfying the necessary condi- 
tions and the trace-preserving condition are reported in Table [21 It is easy to check 
that the case of rank(X/j) = 2 which would be possible only for k = or k = 1 
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doesn't satisfy the iff condition (|56|l . Therefore it is possible to have only rank-one 
operators Xk- 

As a specific optimization problem, let's consider the maximization of the fidelity 
averaged over the two outputs 

F = (^{Tr^d^D} + Tr 2 [^(|^)(V|)]}|V) 

= Tr[i(J® |^><VI + MM ®J>HIVXV>I)] 

and for equatorial qubits we can choose \if>) = |+), where |±) = -^(|0)±|1)). Then 
the fidelity rewrites as 

(66) F = Tr[WR^], 

(67) w = \+){+\ m + Kl-X-I ® l+X+l + l+X+l ® l-X-l) ® l+X+l- 

One can see that W is invariant for permutations over the output copies, and, 
by construction, also all vectors in Table |21 have the same symmetry. Due to the 
special form of the fidelity, the optimal map (satisfying ^{Ij^) = 1 is obtained 
for S = {0, 1} with corresponding rank-two operator given by 

J R^=I^ (0) )^ (0) I + I^ (1) XV' (1) I, 

(68) |^ (0) >- =^(|000) + ^|011) + ^|101», 

|^)=^(|111>| + ^|10Q) + ^|010)). 

6.2. Example. Consider the phase-covariant cloning flfil I22| for equatorial qubits 
from 1 to 3 copies. This correspond to G = U(l), with representations — 
e ^|i><i|o andV^ = e^^Li WWh where s = denotes the input qubit and s = 1,2,3 
the output ones. Here Jff = C 2 and = Jf® 3 . We first need to decompose the 



k 


\h) O \hj) 


-1 


|0001) 





|1001), |0101), |0011), |0000) 


1 


|1000), |0100>, |0010), 1 1101) , |1011), |0111) 


2 


|1100), |1010), |0110), |1111) 


3 


|1110) 



Table 3. Orthonormal bases for the supporting spaces J4?k <8> 
C m ' = C" lfc of the fcth equivalence class of irreducible representa- 
tions for 1 to 3 phase-covariant cloning. The orthonormal basis are 
chosen as subsets of an orthonormal basis for the tensor product 

representation ® U%. This is made of one-dimensional representations, with 
characters e lk ^, with k = —1,0, 1,2,3 and multiplicities rn_i = 1, mo = 4, mi = 
6, mi — 4, and m.3 = 1. Orthonormal bases for the supporting spaces ® 
C" lfc = C" lfc of the fcth equivalence class of irreducible representations are reported 
in Table as subset of an orthonormal basis for the tensor product jff <g> . 
Again, since dim(Jf ) = 2, the necessary condition H64|l says that we can have 
only one equivalence class fc with rank(Afe) < 2, or two equivalence classes both 
with rank(Afc) = 1. In Ref. [22J it is shown that the map which optimizes the 
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averaged equatorial fidelity is actually given by the rank-one map for S = {1} with 
corresponding operator Rj( given by 

i?^H^ (1) W (1) l, 

(69) M , 

=4=(|1000) + |0100) + |0010) + |1101) + |1011) + |0111)). 



Notice that, as a consequence of the specific symmetric form of the chosen fidelity 
criterion, the cloning maps of the examples 16.11 and 16.21 are both symmetrical, 
namely the output Hilbert space is indeed restricted to the symmetric tensor space 
(Jf?® n ) + . Clearly, with the same method also nonsymmetric types of cloning can 
be analyzed well. 

6.3. Example. Consider a generic covariant QO with J(f ~ J^, V g = U g and 
G = SU(d), where d = dim(Jff). In this case the representation U g (8 U* has 
two irreducible components: one which is one-dimensional, corresponding to the 
invariant vector |7) £ Jt?® 2 , and one on the orthogonal complement, and the two 
components will be denoted by k — and k = 1, respectively. Since both the 
irreducible components of the representation have unit multiplicity, the operator 
R = X^X must have X — J2kes °kPk, Ck £ C, Pk denoting the orthogonal projector 
on the invariant space of the irreducible component k, and the necessary condition 
(|64fl is trivially satisfied. On the other hand, one can see that the iff condition 
(|56fl is satisfied for the irreducible representations S = {0} and S = {1}, whereas 
for the reducible one S = {0,1} the map 2?{0) = Tr,^[X^OX] is never injective 
on A' G n B(RngpQ) (one has Ti^[X^OX] = ±[|c | 2 a + (d 2 - l)\ Cl \ 2 ai ]I^ for 
O = ao-Po + fti-Pt) a o^O'i <= C). Therefore, the only trace-preserving optimal maps 
are those corresponding to the operators R = and R = ■^r^i{I® 2 — h\I){I\)i 
corresponding to the trivial map ^ = <# and to the so-called isotropic depolar- 
izing channel .-#(0) = Tr[0]/j^ — P- Finally, notice that in the present 
example the optimal covariant maps are compatible only with (multiple of) the 
trace-preserving condition, since both partial traces Trj^Pfc] are proportional to 
the identity. 

6.4. Example. We consider now the same problem as in the previous example, 
but now with V g = U*. In this case we need to consider the positive operators R 
which are invariant under U* ®U*. It will be easier to consider the representation 
U g (g)U g and then take the complex conjugate of R at the end. Now we have again two 
irreducible inequivalent components, say k = ± with invariant spaces (^® 2 )±, the 
symmetric and the antisymmetric spaces. As in the previous example, the general 
form of R = X^X is X = J2 ke s c k p k, c k € C, and P ± = ± E), where E is 
the swap operator on the tensor product. However, the map &{0) = Trj(f[X'OX] 
is iniective on A' G n £>(Rng(X)) only for representations with a single irreducible 
component. One can see that Tr^[P±] = i(rf± l)Ijt°, and only trace-preserving 
(or multiplying by a constant) QO's are compatible with the present covariance. 
In conclusion, the only extremal covariant operators are R± = (d± 1) -1 (/® 2 ±E), 
corresponding to the channels J%±(0) = (d± l) _1 [Tr(0)/,#' ± O t \. The map 

is the optimal transposition map of Ref . j2EJ ■ 
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